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CN ■ Abstract 
X> ■ 

<]} , We consider the recursion operators with nonlocal terms of special form for evolution systems in (1 + 1) 

pLn ' dimensions, and extend them to well-defined operators on the space of nonlocal symmetries associated with 

, the so-called universal Abelian coverings over these systems. The extended recursion operators are shown 

to leave this space invariant. These results apply, in particular, to the recursion operators of the majority of 
known today (1 + l)-dimensional integrable evolution systems. We also present some related results and de- 
' scribe the extension of them and of the above results to (1 + l)-dimensional systems of PDEs transformable 

. into the evolutionary form. Some examples and applications are given. 



1 — 1 1 Introduction 

| The scalar (1 + l)-dimensional evolution equation possessing an infinite-dimensional commutative Lie algebra 
' of time-independent local generalized (Lie-Backlund) symmetries is usually either linearizable or integrable 
q ■ via inverse scattering transform, see e.g. p] - [El], |8|, 0, [L4|, |18| for the survey of known results and |1C] 



CN '. for the generalization to (2+1) dimensions. The existence of such algebra is usually proved by exhibiting the 
recursion operator [|l2| or master symmetry But in order to possess the latter the equation in question 
• must have higher order time- dependent symmetries, which usually turn out to be nonlocal. This fact is one 
""^h , of the main reasons for growing interest in the study of the whole algebra of time-dependent symmetries of 
• '~| j evolution equations Moreover, for the evolution equations with time-dependent coefficients it is natural 
C ■ to consider their time-dependent symmetries from the very beginning j|, [7| |i~5|| . 

■ Nowadays there seem to exist two basic approaches to the definition of nonlocal symmetries of PDEs. They 
^ | are presented in the papers of Fuchssteiner |3] and the book of Blaszak jlj, and references therein, and in the 

$h ' works of Vinogradov et al., see e.g. || 16, |l7]] and references therein. Our approach is a combination of both. 

' 



In the present paper we shall consider the objects called by Vinogradov et al. the shadows of nonlocal 
symmetries associated with the so-called universal Abelian covering (UAC) over the evolution system (|l]). It 
was shown by Khor'kova || that in the case of UAC any shadow can be lifted to a nonlocal symmetry over 
UAC in the sense of H 16]. Moreover, it is easy to see that the nonlocal symmetries in the sense of [^, 16] 



with zero shadows are in a sense trivial and hardly represent any practical interest. Hence, we lose no essential 
information about the nonlocal symmetries of (|]), when we restrict ourselves to considering just the shadows. 

Since most of authors, see e.g. pi 0, |L [ll| and references therein, call similar objects just "nonlocal sym- 
metries" , in this paper we shall essentially keep this tradition, calling the objects of our study "nonlocal UAC 
symmetries" . This does not lead to any confusion, because nonlocal symmetries in the sense of @, Ul will not 
appear in this paper. The precise definitions and further comments are given in Section ^ below. Note that 
our interest in nonlocal UAC symmetries was inspired by the results from |^, [Tq| , where it was shown that 
in certain cases it is possible to extend the action of the recursion operator so that its application to these 
symmetries produces symmetries of the same kind. 



*The research was supported by Grants CEZ:J10/98:192400002 and VS 96003 "Global Analysis" of the Czech Ministry of 
Education, Youth and Sports and by the Grant No. 201/00/0724 of the Czech Grant Agency. 
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In Section [|| we prove two main results of the present paper. The first one states that a large class of 
recursion operators for systems (|l]) can be extended to well-defined operators, namely to recursion operators 
in the sense of Guthrie ||, acting on the space of nonlocal UAC symmetries and leaving this space invariant. 
The second result states that under certain conditions, which are satisfied for the majority of known ex- 
amples, the repeated application of (extended) recursion operator to weakly nonlocal UAC symmetries (i.e., 
to the symmetries that depend only on the nonlocal variables of the first level, that is, on the integrals of 
local conserved densities, and are linear in these variables) again yields weakly nonlocal UAC symmetries. 
An important consequence of the latter result is that the hereditary algebras (see e.g. ]!]] for definition) of 
time-dependent symmetries for evolution systems (|]) are usually contained in the set of weakly nonlocal UAC 
symmetries. Note that our results are in a sense complementary to Wang's [18] sufficient conditions of locality 
of time-independent symmetries, obtained with usage of recursion operator. 

In Section |I] we explain how the results of Section || can be transferred from (1+1) -dimensional evolution 
systems to arbitrary (1 + l)-dimensional systems of PDEs transformable into the evolutionary form by the 
appropriate change of variables. This is illustrated by the example of sine-Gordon equation. Finally, in Section || 
we discuss further applications of our results and consider the example of Harry Dym equation. 

In Appendices A and B we present two useful technical results: the characterization of kernel of the 
total x-derivative in the space of nonlocal UAC functions and the proof of the fact that the set of nonlocal 
UAC symmetries for (|l]) is a Lie algebra under the so-called Lie bracket — a natural commutator for nonlocal 
symmetries (cf. e.g. [||]). Note that if the hereditary algebra of time-dependent symmetries for (||) is generated 
by (a finite or infinite number of) nonlocal UAC symmetries, then by virtue of the latter result all elements 
of this algebra are nonlocal UAC symmetries. 



2 Basic definitions and structures 

Let us consider a (1 + l)-dimensional system of evolution equations 

du/dt = F(x,t,u,...,u n ) (1) 

for the s-component vector function u = (u 1 , . . . , u s ) T . Here Uj = d^u/dx^ , uo = u; F = (F 1 , . . . , F S ) T ; T 
denotes the matrix transposition. As in |0, |14|| , we make a blanket assumption that all functions below (F, 
symmetries G, etc.) are locally analytic functions of their arguments. This allows, in particular, to avoid the 
pathologies caused by the existence of divisors of zero in the ring of C°°-smooth functions. 



2.1 Universal Abelian covering 

Following |], 16|, describe the construction of universal Abelian covering over a system of PDEs for the 
particular case of system ([l]). 

A function f(x, t, u, uj, . . .) is called local (cf. [||, ||, [ji|), if a) / depends only on a finite number of variables 
Ufc and b) / is a locally analytic function of its arguments. 

The operators of total x- and i-derivatives on the space of local functions are defined as (cf. e.g. 12, Ch. V]) 

Df ] = d/dx + £ E uh x d/dul, 



1=1 i=0 

d{ 0) = d/dt+ £ ziD^yiF^d/du 1 ,. 

1=1 i=0 



Let {Dl°\pa^) = D^'\aa' > ) \ ct £ X{\ be a basis for the space CL^ of nontrivial local conservation laws for 
(|l]) considered modulo trivial ones. Recall that locality means that pa and are local functions and non- 
triviality means that p^ ImD^ (i.e., p^ cannot be represented as a total x-derivative of a local function). 

We introduce || [l6|] nonlocal variables u>a^ of the first level as "integrals" of p^ ■ Namely, we define them 
for all a£li ELS £1 solution of the system of PDEs 

dJ^/dx = p { S\ , . 

dJV/dt = a£\ U 
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It is clear that ui^ are nothing but the potentials for the conserved currents (pa 1 , a a l 

Now let i 
the formulae 



.(0) joh 



Now let us extend the action of operators Ar°^ and to the functions that depend on w^' by means of 



aeli 



.(0); 



X 1 ) 



(i) 



A (0) + E ^ra/ao;, 



,(0); 



(1) 



and consider a basis {d\ (pa } ) 



means that pa $ ImDi 1 ') local conservation laws for the system of equations (||), (H). The densities p a ^ and 



a £2^} in the set CLp of all nontrivial (nontriviality now 

(i) 



fluxes cri may depend not only on x, t, u, Ui, U2, . . ., but on as well, and any given pa' or <Ja' depends only 
on a finite number of u r and of nonlocal variables u„ . Here denotes the totality of variables tc4 f° r « £li. 
We further define the nonlocal variables of second level by means of the relations 

duja^ /dt = aa\ a G X 2 , 

extend the action of A and to the functions that may depend on u a 2 \ and so on. 

Iterating this procedure infinite number of times, we obtain an infinite-dimensional covering hi over ([[]), 
which is called universal Abelian covering (UAC) |], 16]. More precisely, the covering constructed in this way is 
just a representative of the class of equivalent coverings, and the authors of || 16] identify UAC with this class. 



(i) 

Thus, U involves the infinite set of nonlocal variables uj)£ defined by the relations 



dJ^/dt 



,a€lj,j G N, 



Pa 

o-b-V^eljJ €N. 



(*)/■„(*)> 



Here %k+i> fc > 1, is a set of indices such that the conservation laws D± (p a 

(*) 



basis in the set CL^, of all nontrivial local conservation laws of the form A* \p) = D x (a) for (jl|) and (||), (||) 

(k) 

with j < k. The locality means that the densities p and fluxes a of conservation laws from CL]/ depend only on 
x, t, u, ui, . . . and , but not on {3^ m ' with m > k, and any given density or flux depends only on a fi- 

nite number of u r and of nonlocal variables Wa • The nontriviality of a conservation law (p) = (a) from 



(3) 
(4) 

for a G Ik+i form a 



CLp^ means that p cannot be represented in the form D^ ! (f) for some / = f(x,t,uj 



1(1) 



,d;^,u, ui,. 



We employed here the notation 



A 



(fc) 



A 



(fe-i) 



+ E p. 



(fc-1). 



a/a^i fc) , A (fc) = A (fc_1) + E <&~ x) d/dJ£\ 



and the notation uj^ for the totality of variables to a k \ a € Xjt. We shall also denote by cJ the totality of 

(7) 

variables u>a for all j and a. 

The operators of total derivatives on the space of functions of x, t, u, u, Ui, . . . are 



D = D, 



d [ x ] + E E 



A (0) + E E 0-, 



j=l aeJj 

(i-i) 



d/dJJl 



j=l a£jj 



The relations A^(pa^) 



the latter and of the equality [Ac , A 



D<£\<£>) 
(0) n (o)- 



imply the compatibility of 
= are the relations 

[D x k) ,Dl k) ] = 0,k = l,2,..., 



and (01). In turn, the consequence of 



[Ac, A 



0. 



We shall say (cf. §§, [k| [Tl|) that a function / = f(x, t, uj, u, ui, . . .) is a nonlocal UAC function, if a) / de- 
pends only on a finite number of variables Wq' and and b) / is a locally analytic function of its arguments. 
We shall call / a nonlocal UAC function of level k, if / is a nonlocal UAC function independent of uia for j > k. 



3 



Since the kernel of D in the space of nonlocal UAC functions is exhausted by functions of t (see Appendix 
A for the proof) , it is easy to verify that our definition (cf . || ) of nontriviality of a conservation law is in fact 
equivalent to the standard one |l2|, |l6| . 

Let us stress that p£ ^ are defined up to the addition of the terms from ImD^ 1 \ and a a ^ are defined 

(k—l) 

up to the addition of the terms from ker Dx , i.e., they should be considered as equivalence classes modulo 
ImD x k ^ and hev D x k , respectively. This means, in particular, that the nonlocal variables are defined 
up to the addition of arbitrary nonlocal UAC functions of level k — l. 

Making different choices of representatives in these equivalence classes yields different coverings over (|l|) , 
but these coverings are equivalent in the sense of definition from [jO], Ch. 6], and in the sequel we shall assume 
that we deal with a fixed representative of the respective equivalence class, because constantly operating with 
the whole class in the explicit computations is extremely inconvenient. However, the results obtained below 
are obviously independent of this choice, and thus hold true for the whole class of equivalent coverings, which 
result from the above construction. 



2.2 Nonlocal UAC symmetries 



We shall call (cf. 11, [T^, ^6[) an s-component nonlocal UAC vector function G a nonlocal UAC symmetry 
of (|l|), if the evolution system du/dr = G is compatible with (|l|), i.e., d 2 u/dtdr = d 2 u/drdt, where the 
derivatives with respect to t and r are computed with usage of ([!]), (||) and du/dr = G, respectively. We 
shall denote the set of all nonlocal UAC symmetries of (jlj) by NSf(U). If dG/duj = 0, then G is called (local) 
generalized (or higher local, or just local) symmetry of (|l|), see e.g. [||, ffO), [l^]. With G being a nonlocal UAC 
vector function, the compatibility condition for ([l|) and du/dr = G takes the form 

A(G) = F'[G], (5) 

where F' = £JL d¥/dUiD\ 

Let us mention that nonlocal UAC symmetries are nothing but a particular case of general nonlocal sym- 
metries, considered e.g. in pi 0. Indeed, the determining equations for the latter, given in (l], §], are nothing 
but the compatibility conditions for ([!]) and u T = G, and these conditions reduce to (||), if G is a nonlocal 
UAC vector function. 

The set ~NSf(U) of nonlocal UAC symmetries is interesting and important. In particular, our results imply 
that for nearly all known examples the action of the recursion operator is well defined on NSf(U) and leaves it 
invariant. Hence, ~NSf{M) contains all elements of the hereditary algebra (see ||], ||, [| for its precise definition) 
of time-dependent symmetries for (H), if they are generated (cf. e.g. Jlj) by means of the repeated application 
of the recursion operator to the scaling symmetry of ([!]) and to a time-independent local generalized symmetry 
of © (e.g. to F, if dF/dt = 0). Note that NS F (U) is a Lie algebra with respect to the so-called Lie bracket 
(see Appendix B below for the proof). Therefore, if the hereditary algebra for ([!]) is generated by (a finite or 
infinite number of) nonlocal UAC symmetries, then all its elements are nonlocal UAC symmetries. Note that 
there also exist (see e.g. M, 15]) integrable systems (|]) that possess only a finite number of local generalized 



symmetries, but have infinite hierarchies of nonlocal symmetries, and these nonlocal symmetries turn out to 
be nonlocal UAC ones. 



To avoid possible confusion, let us stress that the definition of nonlocal symmetries used in [y, 16] is 



different from the ours, but by Theorem 3.1 from |q] we always can recover from G G NSi?(W) the nonlocal 



symmetry in the sense of [o, 16 [. 



3 On the action of recursion operators 

In this section we prove for the case of systems (|l|) and recursion operators of the form ([?]) the conjecture 
of Khor'kova || stating that for any (1 + l)-dimensional system of PDEs its recursion operator 91 can be 
extended to a well-defined operator 9\ on the space NSf(U) of nonlocal UAC symmetries and leaves this space 
invariant. 
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In complete analogy with the case of local functions, see e.g. [12, Ch. V, §5.3], we call the operator 

i . <i . 

= ^2(—Dy o of a formal adjoint of <B = ^ biD 1 . Here hi are some r x r matrix-valued nonlocal UAC 

i=0 i=0 
functions, T denotes the matrix transposition and o stands for the composition of operators. 

An s-component nonlocal UAC vector function 7 is called a cosymmetry jj], [l8j of (||), if it satisfies the 
equation 

A(7) + (F') f [7]=0. (6) 



Nearly all known today recursion operators for systems (||) have the form (cf. [jl8|| ) 

k p 
i=0 3=1 

where Oj are s x s matrix-valued local functions, and Gj and -jj are local symmetries and cosymmetries of 
(||), respectively. 

Note that the action of recursion operators of the form ([?]) is initially defined only on local generalized 
symmetries (see e.g. [|l^]), but the formula (0) together with the subsequent definition of D -1 enable us to con- 
struct an extension 9i of 91 to the space NS.p(W). However, it is not clear a priori whether £K is a well-defined 
operator on ~NSf{M) and whether SH(G) is a nonlocal UAC symmetry of ([[]), provided so is G. 

Let us mention the following result of Wang |[^]. Suppose that 9^ (Q), with Gj and 7,- being arbitrary 
s-component time-independent local vector functions, is a recursion operator for (Q), both 91 and system ([lj) 
are homogeneous with respect to a scaling of x,t and u, and dF/dt = and dJi/dt = 0. Then under minor 
restrictions on F and 9^ the functions Gj and -fj indeed are symmetries and cosymmetries for (|l]), respectively. 
Obviously, this result remains valid when Oj, Gj and 7,- are nonlocal UAC functions. 

In order to proceed, we should agree how to interpret the action of D^ 1 . For our purposes it suffices to 
adopt the following definition, which is a particular case of the general construction of Guthrie @: 

Definition 1 Let P be a nonlocal UAC function such that Dt(P) = D(Q) for some nonlocal UAC function 
Q (i.e., P is a conserved density). 

Then we shall understand under R = D^ 1 (P) a solution of the system 

D(R) = P, ^ 
D t (R) = Q. W 

The existence of nonlocal UAC solution R of @, provided D t (P) = D(Q), is one of the fundamental 
properties of universal Abelian covering, proved in ||. Note that according to the above definition we have 
Ua = P > l {pa as it would be natural to expect. As we show in Appendix A, the kernel of D in the space 
of nonlocal UAC functions is exhausted by functions of t, whence it is immediate that the solution R of (||) 
for given P and Q is unique up to the addition of an arbitrary constant. So, the "integral" D^ 1 (P), defined 
above, is not an integral of P itself, but rather of a conservation law Dt(P) = D(Q). It is also clear that the 
integrals D~ 1 {P) evaluated for different Q in general will differ by a function of t. 

In order to prove that we can extend the action of 9\ to any nonlocal UAC symmetry H of (|l]), we have 
to show that for any H G NSp(U) we have Dt(~fjH) = D(Cj) f° r some nonlocal UAC functions Q. Indeed, 
then the integrals D~ 1 (7 J H), interpreted in the sense of the above definition with Q = Q, are nonlocal UAC 
functions defined up to the addition of arbitrary constants. Note that the expressions like ab stand here and 
below for the scalar product of s-component vectors a and b. 

By © and (|) we have A(H) = F'[H] and A(7j) = -(F'^fy], whence 

A(7,H) = -(F') t [7 J ]H + 7j F / [H]. (9) 



There is an obvious generalization (cf. e.g. [12, Ch. V, §5.3]) of the well-known Lagrange identity from 
theory of ordinary differential equations, namely 

Mg) - &(f)g = D( V ), 7]=J2 Ei-DyibJf^D^-^g), (10) 

i=lj=0 
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<? 

valid for any differential operator 23 = Yl ^D 1 and for any r-component nonlocal UAC vector functions / 

i=o 

and g, provided b% are r x r matrix-valued nonlocal UAC functions. 

Using ([n]) for 23 = F', / = 7j, g = H, we conclude that £^( / y J -H) indeed can be represented in the form 
D(Q) for the nonlocal UAC function 

n t— 1 

= EE(- I) ) m (( 5F /^)^)^ 1 ( H )' ( n ) 

i=l m=0 

Hence, for P = 7jH we can always make a 'canonical' choice Q = Q while computing D~' 1 (P) according 
to Definition |[ With this choice and the above definition of D^ 1 , the recursion operator 21 (f^) is, in essence, 
replaced by a new operator 21, which is easily seen to be a recursion operator in the sense of Guthrie pj]. 

Many important properties and definitions (for instance, that of hereditarity) can be readily transferred 
from 21 to 21. The operator 21 is free Jl of the pathologies caused by nai've definition of -D -1 , cf. [13] for an 



alternative way of overcoming these difficulties. In particular, it is easy to see that 21 always maps nonlocal 



UAC symmetries to nonlocal UAC symmetries (cf. ||]), because 21 satisfies the equation [D t — F',21] = [12]. 
Thus, we have proved 

Proposition 1 Any recursion operator 21 for ^), with a% being s x s matrix-valued nonlocal UAC func- 
tions, and Gj and~/j being nonlocal UAC symmetries and cosymmetries for ([ip, respectively, can be extended 
to a well-defined operator 21 that acts on the whole space NS^(W) of nonlocal UAC symmetries for ffi) and 
leaves this space invariant. 

Now let us consider what happens if cij, Gj,~fj are local and we apply 21 to a local generalized symmetry 
H of (|l]). The above reasoning indeed holds true. Moreover, we see that 7jH are local conserved densities for 
(H) and the respective Q are local as well. Hence, the application of 21 to local generalized symmetries of (||) 
yields nonlocal UAC symmetries of the form (cf. || and fl6| , Ch. 6]) 

H = H +]TH aW «, (12) 

ael H 

where Ho and H Q are s-component local vector functions, and Ih is a finite subset of X\. 

We shall call an r-component nonlocal UAC vector function weakly nonlocal UAC vector function, if it can 
be represented in the form (12) with Hq and H Q being r-component local vector functions. We shall denote by 



WNLSi?(ZY) the set of all weakly nonlocal UAC symmetries for (|l]). For the majority of integrable systems (||) 
their master symmetries are s-component weakly nonlocal UAC vector functions. 
Note that if H is a local generalized symmetry of (|l|) and 

7i H = (13) 

where £j are local functions, then by the above Z) _1 (7jH) can differ from £j only by a function of t. Hence, 
D^^jH) is a local function, and thus 21(H) is a local generalized symmetry for (|l|). In other words, the appli- 
cation of 21 to local generalized symmetries of (|l]) satisfying ( |l~3|) again yields local generalized symmetries of 
(|). Below we shall assume (obviously without loss of generality) that Gj in (0) are linearly independent. Then 
it is easy to see that the conditions (|l3l) for j = 1, . . . ,p are equivalent to the requirement that 21(H) is a local 
generalized symmetry, provided so is H. Let us mention that Theorems 6-8 and 6-9 of Wang pi] provide an 
easy way to verify the conditions ([D^) for large families of time-independent local generalized symmetries of (Q). 

Since F and the coefficients of D+ are independent of u>a for all a and j, it is immediate that for any 
nonlocal UAC symmetry G of level one for (|l]) the quantities dG/dtOa^ satisfy the determining equation (p) 
and hence also are nonlocal UAC symmetries of level one for (|l]). In particular, for any H of the form fll2 ) 
the quantities H a = dH/duJa^ are in fact local generalized symmetries of ([l]). 

Using this result, let us show that H = 21(H) € WNLS F (W) for any H of the form |[|), provided 2t(H Q ) 
are local generalized symmetries of (|l]) (or, equivalently, 7jH Q = D(£j >a ) for all j = 1, . . . ,p and all a £ Ih, 
where £j >a are local functions). 



6 



As D^oj^) = D i - l (p { a ) ) are local functions for i > 1, it is clear that <K(H) G WNLS F (^), if there exist 
scalar weakly nonlocal UAC functions Rj such that 

D(R J )= lj U, (14) 
D t (R j ) = C j , (15) 

where C,- are given by (|Tl|). Indeed, then fR(H) is a weakly nonlocal UAC vector function, and by Proposition |] 
91(H) e NS^(W), hence «K(H) G WNLS F (W). 

Let Rj = Rj- £ Zj^J^- Using (g), (|l|), we obtain 

D(Rj)= lj U - £ e jiQ ^ 0) =Vi, 
A(^) = Ci- £ A(£i,aO;«) = Xi . 

It is clear that V'j are local functions and that Dt(ipj) = D(\j)- If we show that Xj are local as well, then 
D~ 1 (ipj) obviously are linear combinations of (modulo local functions), so Rj are weakly nonlocal UAC 
functions, and the result follows. 

As Xj may depend on the nonlocal variables of the first level tc4 at most, we only have to check that 
dxj/duJa^ = 0, i.e., d(j/duia^ = D t (£,j : a). Obviously, the only nonlocal terms in Q are 

E ( EE(- Z} ) m (( aF /^) T ^) Z)4 ~ m " 1 ( H «) w «' 

ael H i=l ™=0 ' 

Hence, in order to prove our result it remains to show that 

n i— 1 

Ci,a = E E(- J D) m ((9F/9u i ) T 7i ) J D i - m - 1 (Ha) = Afe,a). 
i=l m=0 

Comparing this equality with (11) and bearing in mind that H a are local generalized symmetries of (ffl), 
we see that 

= A(7,H a ) = = f(A(^,a)). 

Using Proposition || from Appendix A, we find 

Cj,a = Cj,a{t) Dt{£,j,ot) j 

where Cj ja (i) is arbitrary function of i. 

But it is clear that the function ^ iQ , determined from the relation D(£j )Q ) = 7jH a , is defined only up to 

the addition of arbitrary element of ker D, i.e., an arbitrary function bj tOC (t) of t. Hence, replacing £j a by — 

t 

J Cj ia (T)dr, we can assume without loss of generality that Cj )(X (t) = 0, and thus Cj,a = D t (£,j t0l ), as required. 

to 

Thus, we have proved the following result, generalizing Proposition 4.1 from Q: 

Proposition 2 Let ^) have a recursion operator $H (Q), where ai,Gj,~fj are local. Then for any H G 

WNLSf(^) the quantity SK(H) is well defined and <H(H) € WNLS F (^), provided ^(dH/du^) are well- 
defined local generalized symmetries of ^) for all a G T\ . 

If H = £H(G) for some local generalized symmetry G, then H Q in (12) are linear combinations of the 
symmetries Gj that enter into 9\. We can easily see that the coefficients H Q = d"H./doJa at oj^ in the 
representation ( |i~2| ) for H = 9v(G) are linear combinations of 5H(H a ) and of Gj, and hence are in fact linear 
combinations of Gj and $R(Gj) only. Thus, H G WNLS_f(^0, provided !SH(Gj) are well-defined local general- 
ized symmetries of (|l|). This is equivalent (cf. above) to the requirement that 7j!SK(Gi) = D(£ij), where £y 
are local functions, for all i,j = l,...,p. 
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Iterating this reasoning, we conclude that JH(Q) can be represented in the form (|7|) for any nonlocal UAC 
symmetry Q obtained by the repeated application of the recursion operator $H to local generalized symmetries, 

provided T^JH^Gj) = Dfaj^), where Ci,j,d are local functions, for all i,j = 1, . . . ,p and all d = 0, 1, 2, 3, 

These conditions are equivalent (see above) to the requirement that 9\ d (Gi) are well-defined local generalized 
symmetries of (||) for all i = 1, . . . ,p and d = 1, 2, 3, . . .. 

In particular, if these conditions are satisfied, then for any (time-dependent) local generalized symmetry 
G of (|) we have VV(G) € WNLS_f(W) for all j = 1,2,.... Hence, if the hereditary algebra (see e.g. @ for its 
definition) of time-dependent symmetries for (0) is generated by the repeated application of the extension 9\ 
of a recursion operator 9\ (0) to some local generalized symmetries of (|l]), and dj, Gj and ~fj are local and 
satisfy the above conditions, then all elements of this algebra belong to WNLSir(W). 



4 Generalization to non-evolution systems 

The above results can be applied to any (l + l)-dimensional systems of PDEs transformable into the evolution- 
ary form (H) by the appropriate change of variables. This set includes, in particular, all systems transformable 
into Cauchy-Kovalevskaya form 

^ = ^(x, t, u\ . . . , u\ . . . , d a+ ?u J /dt a dxP, ...), I=l,...,q, (16) 
where 3>/ may depend only on x, t, u , . . . , u q and 



\d a+p ujldt a dxP\a <rj — 1,(3 <k\ , J = 1, . . . ,q. 



The system (|T^) can be further transformed into an evolution system of the form (0) by introducing new 
dependent variables v l a = d'^u 1 jdt a for a = 1, ... ,77 — 1. Indeed, combining the variables u 1 , . . . , u q and 
into a single vector v, we see that (jig) together with the equations 



/8t a =v T a , a = l,...,n-l,I = l,...,q, 
forms the evolution system of exactly the same form as (|l|): 

dw/dt = K(x, t, v, dv/dx, d 2 v/dx 2 , d k \/dx k ). (17) 



The class of systems of PDEs transformable into the form fllq ) and hence into the evolutionary form (17) 
is very large. In particular, it includes |l^, Ch. 2] all analytic locally solvable (1 + l)-dimensional systems of 
PDEs possessing at least one noncharacteristic direction. The majority of known examples of non-evolutionary 
integrable (1 + l)-dimensional systems are indeed transformable into the form (f[6|) by the (appropriate mod- 
ification of) above change of variables. 

Hence, Khor'kova's conjecture stating that the (extended) recursion operators are well defined on nonlocal 
UAC symmetries holds true not only for the evolution systems (II]) with the recursion operators (0) , but also 
for any systems transformable into Cauchy-Kovalevskaya form ( |l6|) and then into (|l7|) , provided the recursion 
operator for transformed system (|l7]) has the form (^). Indeed, making the inverse change of variables we can 
readily see that the (extended) recursion operator of original system is also well defined on its nonlocal UAC 
symmetries. 

For instance, it is well known that the sine-Gordon equation u^ v = sinu can be transformed into uu — u xx = 
sinu by setting x = ^ — 77, t = Then, introducing a new dependent variable v = ut, we obtain the evolution 
system (see e.g. [^] and references therein), equivalent to the SG equation: 

ut = v, v t = u xx + sinu. (18) 

The recursion operator 9\ for (|l^) is (see e.g. [18|) of the form (0), so by Proposition |l| the action of 



9\ on nonlocal UAC symmetries of ( |18D is well defined and leaves the space of these symmetries invariant. 
Returning to the original variables, we conclude that the same is true for the recursion operator (rewritten 
as a recursion operator in the sense of Guthrie ||]) of SG equation, so we recover the result of Khor'kova ||, 
initially obtained by straightforward computation. 
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5 Applications 



Consider, for instance, the well known integrable Harry Dym equation ut = v?u^. Its recursion operator (see 
e.g. |l8|) 9\ = u 2 D 2 — uu\D + UU2 + ^u^D^ 1 o u~ 2 indeed has the form ([?]), and u~ 2 is a cosymmetry and 



u u% is an (obvious) symmetry for this equation. Using (11), we find that D t (u G) = D(D (uG) — 3u\D(G)) 



for any nonlocal UAC symmetry G of HD equation, and hence D~ l (u~ 2 G) is a nonlocal UAC function. Thus, 
by Proposition [l] the extension 9t of the above 9\ is well defined on the space of nonlocal UAC symmetries of 
HD equation and leaves this space invariant. 

It is possible to give a lot of other examples where Proposition [I] ensures that the (extended) recursion 
operators are well defined on the space of nonlocal UAC symmetries and leave it invariant. This fact often 
allows to draw a number of useful conclusions. For instance, provided the hereditary algebra |l], ||] of time- 
dependent symmetries is generated by the repeated application of the extension 91 of a recursion operator 
JR of the form (fl) to some local generalized symmetries, all elements of this algebra are nothing but nonlo- 
cal UAC symmetries. Moreover, if 91 fc (Gj) are well-defined local generalized symmetries for all k G N and 
j = 1, . . . ,p, then by Proposition [2] all elements of this algebra are weakly nonlocal UAC symmetries, that is, 
they depend only on the nonlocal variables u>a \ i.e., on the "integrals" of nontrivial local conserved densi- 
ties for (|l]), and are linear in u„ . The hereditary algebra contains time-dependent symmetries of arbitrarily 
high order, hence their directional derivatives are formal symmetries of arbitrarily high order (see e.g. (l(J 
for definition of formal symmetry) for (|l]). Thus, the evolution systems (|l|) possessing the hereditary algebra 
have time-dependent formal symmetries of arbitrarily high (and hence of infinite) order, and the coefficients 
of these formal symmetries are usually nonlocal (more precisely, weakly nonlocal) UAC functions. 

To conclude, let us mention that it would be very interesting to generalize the results of the present paper 
to the evolution equations with constraints, introduced in [jil] ]. 

Acknowledgements 

I am pleased to express deep gratitude to Dr. M. Marvan for the numerous and highly stimulating discussions 
on the subject of this paper. 

Appendix A: On the structure of ker D 

The aim of this appendix is to prove that the kernel of operator D in the space NL^(W) of nonlocal UAC 



functions consists solely of functions of t, cf. [11] 



Let Aq be the algebra of all scalar local functions under the standard mulplication. 

We shall call an algebra A of scalar nonlocal UAC functions (under standard multiplication) admissible, 
if it has the following properties: 

• for any locally analytic function h(yi, . . . , y p ) and any ctj G A we have 
h(ai, • • • , a p ) G A; 

• A is closed under the action of D and D t ; 

• A is obtained from the algebra Aq by means of a finite sequence of extensions. 

The third property means that there exists a finite chain of admissible algebras Ao, A\, A2, ■ ■ ■ , A m = A 
such that Aj is generated by the elements of Aj-i and just one new nonlocal variable Q = D~ l (r]j), where 
r]j G Aj-i is such that r\j \m.D\j^ ] x and D t {r\j) G ImD| i 4.._ 1 . 

Consider a nonlocal UAC function /. It may depend only on a finite number m of variables Wq , and it is 
easy to see that there exists a minimal (i.e., obtained from Ao by means of the minimal possible number of 
extensions) admissible algebra K, of scalar nonlocal UAC functions, which contains /. 

It is clear that in order to prove that / G ker D implies that / depends on t only it suffices to prove that 
ker D |k; consists solely of functions of t. 

In order to proceed, we shall need the following 
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Lemma 1 Let A be an admissible algebra, ker D\ A consist solely of functions oft, and A be the extension 
of A obtained by adding the nonlocal variable £ = D^ 1 ^), where 7 € A is such that 7 ImDj^ and 

A(7)eIm£»U- 

Then A is admissible and kerZ)|^ also consists solely of functions oft. 



Remark 1. The conditions 7 ImL>|^4 and -^4(7) E Im.D|_4 imply that 7 is a linear combination of p, 
(modulo the terms from ImD\ A ). 



(?) 
a 



Remark 2. This lemma is a natural generalization of Proposition 1.1 from [11] to the case of time- 
dependent nonlocal UAC functions, and its proof relies on the same ideas. 

Proof of the lemma. The admissibility of A is obvious from the above, so it remains to describe ker D\^. 

By definition, the elements of A may depend only on a finite number of nonlocal variables Cii - ■■■> Cm- 
Let Bo C A be the algebra of all locally analytic functions of x, t, u, Ui, . . . , u p , £1, . . . , ( m , where p is the 

minimal number such that Bo contains 7 and -D(Ci) f° r i = 1, ■ ■ ■ ,m. It is straightforward to check that such 

p does exist. 

Consider the following chain of subalgebras of A: 

B j+1 = {he Bj I D\ A (h) = gi,ge Bj}, j = 0, 1, 2, . . . . 

Any locally analytic function of elements of Bj obviously belongs to Bj. 

As Bo is generated by s(p + 1) + m + 2 elements x, t, ui, . . . , ui tP , (±, . . . , Q m , where / = 1, . . . , s, we 
conclude that Bj + \ = Bj for j > s(p + 1) + m + 2 (cf. fn}). Indeed, by construction Bj+i C Bj, and hence 
the functional dimension dj + \ of Bj+i does not exceed that of Bj. If these dimensions coincide, then we have 
Bj + \ = Bj, and otherwise dj + \ < dj — 1. Since do = s{p + 1) + m + 2, it is clear that d = (i s ( p+1 ) +m+2 < 1 
provided Bj + \ 7^ Bj for j = 0, . . . , s{p + 1) + m+ 1. On the other hand, d > 1, because in any case -B s ( p+1 ) +m+2 
contains the algebra of functions of t, and the result follows. Thus, B s ^ p+ i^ +m+ 2 is the algebra of all locally 
analytic functions of some its elements z\, . . . , Zd, i.e., it is generated by Z\, . . . , z&. 

Let / = f(x,t,u, . . . ,Uq,Ci, . . . ,Cm ; C) £ kerL>|^, df /dC, 7^ 0. Differentiating the equality D(f) = with 
respect to Uj, j > p, we readily obtain df /duj = for j > p. Therefore, f £ Bo for any fixed value of C- 

We have 

D(f) = D\ A (f) +1 df/d( = 0. (19) 



But ([L9|) implies that / E B\ for any fixed value of Q. Then, again by virtue of (|19D, we have / G B2 for any 
fixed C, and so on. 

Thus, / G £ s (p +1 ) +m+2 for any fixed value of C and D\ A (f) / 0. Hence, the operator -f% s{p+1)+m+2 is 
nonzero and -D|B s(p+1)+m+2 = 7A, where X is a nonzero vector field on the space of variables Z\, . . . , z^- 

Let w £ S s ( p+1 ) +jn+2 be a solution of equation X(w) = 1. Then D\ A (w) = 7, what contradicts the 
assumption that 7 ImD|_4. The contradiction proves the lemma. □ 

The desired result about kerD|/c readily follows, if we successively apply the above lemma for A = Ao 
and A = Ai, then for A = A\ and A = A2, and so on, until we see that ker D\ Am , A m = fC, consists solely of 
functions of t. 

Our reasoning applies to any nonlocal UAC function /, so we have proved 

Proposition 3 The kernel of the operator D in the space NLj?(W) of nonlocal UAC functions consists solely 
of functions oft. 

From this result it is immediate that the intersection ker D n ker Dt in the space of nonlocal UAC functions 



consists solely of constants, and hence by Proposition 1.4 from 16, Ch. 6,§1] universal Abelian covering over 
(!]) is locally irreducible. 

Let us stress that the above proposition is not valid for the functions that depend on the infinite number 
of variables Uj and lo^ at once. 

Indeed, consider the well-known Burgers equation ut = U2 + uu\, whose only nontrivial local conserved 



density is u (see e.g. [16| for proof). Let ip = ip(x, t, u, u\, . . .) be an arbitrary infinitely differentiable local 
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function. Then it is straightforward to check that the function ^ = ^ n j wnere w = and 

j=0 V- 

T = —(l/u)D, belongs to ker D. It is clear that depends on an infinite number of variables Uj = d 3 u/dx 3 , 
j = 0,1,2,..., provided D(if>) ^ 0. 

Appendix B: Lie algebra structure of NS^(W) 

In this appendix we prove that the set NSp(U) is a Lie algebra with respect to the so-called Lie bracket (see 
e.g. [||), defined as 

[G, H] = H'[G] - G'[H]. (20) 

We employed here the notation /'[H] = (df(x, t, u + eH, ui + eD(H), . . .)/de)\ e= o for the directional derivative 
of any (smooth nonlocal) function / along H, see e.g. [[ij. The bracket ( p0| ) is obviously skew-symmetric. It 
satisfies the Jacobi identity by virtue of properties of the directional derivative, see e.g. [|J and references 
therein. 

00 df ■ 

If / is a local function, then / [H] = ^ — — D l (H) is a well-defined nonlocal UAC function for any 

i=0 oo Qj 

s-component nonlocal UAC vector function H, and f'=Yl ~ — D % is a differential operator, cf. e.g. |p|- (l^1 . 
If / is a nonlocal UAC function, then we have 1=0 % 



= E Z- Dk w + EE ffe^w ( 21 ) 

k=o OUk i=l duX' 



k=0 j=l a£Z. 

Hence, in order to show that /'[H] is a nonlocal UAC function for any H £ NSf(U), we should define u/^[H] 
and show that u'^' [H] are nonlocal UAC functions. 

From now on we assume that the integration constants arising while computing L> _1 according to Defi- 
nition [l] are chosen so that for any constant c we have R = D^ 1 (P) = cR = cD^ 1 (P), where P = cP (and 
Q = cQ). 

Using the interpretation of JJ ] as D' l {p ( i~' 1) ), given above (cf. |, |, 0]), we set (J® [H] = D~ 1 (p'^~ 1 ^ [H]). 
The quantities //^"^[H] can be computed inductively. Indeed, p$ are local functions, so we know the formula 
for p'^[H] (see above), and hence we can evaluate o/^[H]. Next, as p'a 1 involve only uj^ and local variables 
x, t, u, Ui, . . ., we can find p'^ [H], using (|2l|). Then we find p'^ [H], and so on. 

According to our definition of D^ 1 , in order to guarantee that D~ 1 (/9 / ^[H]) are well defined we have to 
show that there exist nonlocal UAC functions such that 

D t (p'^lH]) = D(&). (22) 

As Dtip^) = D(a t i ) ), we have (A(/?i J } ))'[H] = (D{p$))'\H\. Since H G NSj?(W), we readily obtain from (|) 
that (D t (p^)Y[B] = D t (p' { £[H}). It is also easy to see that (£>(fji j) ))'[H] = D(a f{ £[K\). Hence, D t (p' { £[U}) = 
D(a' { J ] [H}), and © holds, if we take <r'^[H] for C,& \ Using this result, we always can make a 'canonical' 
choice Q = cr'^[H] while computing D~ l {P) for P = p'^'[H] according to Definition [j], and thus get rid of 
possible pathologies. 

It is clear from the above that if G,H G ~NSf(U), then G'[H] and H'[G] are nonlocal UAC vector 
functions, and thus [G, H] is a nonlocal UAC vector function. Using (|5|) for G and H, we can easily show that 
[G, H] G NSf(U), so we have proved the following 



Proposition 4 The set NSf(U) is a Lie algebra under the commutator (2C). 
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